JSC-08678 


N/ SA TEC.i .iCAL MEMORANDUM 


NASA TM X- 58126 
January 1974 



AN EFFICIENCY STUDY ON OBTAINING THE MINIMUM WEIGHT 
OF A THERMAL PROTECTION SYSTEM 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
LYNDON B. JOHNSON SPACE CENTER 
HOUSTON, TEXAS 77058 


<NASA-TM-X-581^b) AN EFFICIENCY STUDY ON N74-17644 

OBTAINING THE MINIMUM WEIGHT OF A 
THEFMAL FROTECTION SYSTEM (NASA) 12 p 

HC 14.00 CSCL 20/13 Unclas 

G3/33 3C459 


NASA TM X-58126 


AN EFFICIENCY STUDY ON OBTAINING THE MINIMUM WEIGHT 
OF A THERMAL PROTECTION SYSTEM 


S. D. Williams 

Lockheed Electronics Company, Inc. 
Houston, Texas 

D. M. Curry 

Lyndon B. Johnson Space Center 
Houston, Texas 



AN EFFICIENCY STUDY ON OBTAINING THE MINIMUM WEIGHT 

OF A THERMAL PROTECTION SYSTEM 

By S. D. Williams* and D. M. Curry 
Lyndon B. Johnson Space Center 


SUMMARY 


This report results from an efficiency study conducted to obtain the minimum 
weight of a thermal protection system. Based on previous experience with the 
adaptive creep-pattern search method for the evaluation of thermal protection system 
models, this study sought more efficient methods for reducing computer usage time: 
these methods were a numerical optimization method and two nonlinear least-squares 
methods. Three minimizing techniques were evaluated in an effort to solve the 
minimum -weight problem involving multimaterials and two or more constraints . The 
techniques were the Powell optimization routine, the Powell least-squares routine, 
and the Peckham least-squares routine. Final analysis of data showed that the 
Peckham method of least squares was the most efficient minimizing technique inves- 
tigated , whereas a simple quadratic fit was the most efficient method for solving the 
one material and one constraint problem . 


INTRODUCTION 


The design of a minimum-weight thermal protection system (TPS) for the 
Space Shuttle is desirable to allow the maximum payload capability . This requires 
an analysis involving material selection, arrangement, and thickness subject to 
design constraints of the materials. In a previous study (ref. 1) , the adaptive creep- 
pattern search method was found effective in determining the minimum weight of a 
multilayer insulative TPS with airgaps and two or more constraints . By using this 
analytical approach . the TPS minimum -weight problem is simplified to correlating 
optimum solutions for different TPS configurations . 

Experience with the adaptive creep-pattern search method for the evaluation 
of many TPS models, materials, and trajectories has shown that while the analyst's 
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effort is greatly reduced, tho computer time requirements may be excessive. An 
analysis was therefore conducted to determine a more efficient method for reducing 
the computer usage time and cost . The methods used in this investigation were a 
numerical optimization method developed by Powell (ref. 2) and two nonlinear least- 
squares methods developed by Powell (ref. 3) and Peckham (ref. 4) . As with the 
adaptive creep-pattern search , none of the methods required the formal evaluation of 
derivatives . Finally , so that a meaningful comparison could be made , the same TPS 
models used for the efficiency study in reference 1 were used for the efficiency study 
in this investigation . 

As an aid to the reader , where necessary the original units of measure have 
been converted to the equivalent value in the Systeme International d'Unites (SI) . 
The SI units are written firut, and the original units are written parenthetically 
thereafter . 


ANALYSIS 


The TPS weight optimization problems can be divided into two categories; 
those in which only one material thickness is to be optimized and those in which more 
than one material thickness is to be optimized. The current work showed that for 
the one material and one temperature constraint problem , a simple quadratic fit of 
temperature as a function of thickness is the most efficient technique. Desired tem- 
peratures are calculated for three different trial thicknesses of the material to be 
optimized, and a quadratic curve fit is used to predict thickness as a function of 
temperature. This allows for the direct estimation of a new thickness that will pro- 
duce the desired temperature . The technique is iterative in that , if the temperature 
obtained from the predicted thickness is not within the desired tolerance , one of the 
old thickness values is thrown away and a new thickness is predicted from the 
three remaining values . 

The case for more than one material to be optimized can be further subdivided 
into two classes: when the number of materials to be optimized is greater than the 
number of constraints and when the number of constraints equals or exceeds the 
number of materials to be optimized . The first class requires the use of true opti- 
mization techniques, but, because of TPS design criteria, this class seldom arises 
in practice . The second class is exact or overdetermined , and the solution can be 
made tractable by a least-squares method. In a least-squares approach, no attempt 
is made to achieve a minimum weight. Rather, the problem is to satisfy all con- 
straints; these constraints are temperature constraints imposed on the backwall of 
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the materials in the TPS . In the least- squares formulation , a vector of independent 
variables x = is defined for which a minimum value of S is to be 

determined . 
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with - T* where is the temperature at the backwall of material j and 

T* is the desired temperature for the kth constraint . 


The three minimizing techniques used in this investigation were Powell's VA04A 
optimization routine (ref. 2) , which is a method for finding the minimum of a function 
of several variables without calculating derivatives; Powell's VA02A least-squares 
routine (ref. 3) , which is a method for minimizing a sum of squares of nonlinear 
functions without calculating derivatives; and Peckham's MINSQ least-squares rou- 
tine (ref. 4) , which is also a method for minimizing a sum of squares without calcu- 
lating gradients. 

Three specific TPS configurations were considered for this study: configura- 
tion 1 — Dynaflex/TG- 15 000/ aluminum; configuration 2 — Dynaflex/airgap/ 

TG-15 000/aluminum; and configuration 3 — Dynaflex/TG- 15 000/airgap/aluminum . 
(Dynaflex and TG-15 000 are product names of thermal insulation material.) The 
three configurations , shown schematically in figure 1 , are the same TPS models used 
for the efficiency study in reference 1 . 

The control variables for optimization were the thicknesses of the Dynaflex 
and the TG-15 000. A thin metallic skin was assumed to form the aerodynamic sur- 
face to which the Dynaflex was attached; however, this thin skin was ignored in the 
thermal model. The thickness of the aluminum was held to a constant value of 
0.18 centimeter (0.07 inch) . The maximum temperature at the backwall of the Dyna- 
flex was constrained to 922 ± 0.28 K (1660° ± 0.5° R) , whereas the backwall of the 
aluminum was constrained to 366 ± 0.28 K (660° ± 0.5° R) . No temperature constraint 
was imposed on the backwall of the TG-15 000. The emissivity for the airgap was 
held at 0.2 on both the front and back of the materials. An emissivity of 0.8 was used 
at the face of the Dynaflex for radiation to space . 

Because the purpose of this effort was to evaluate minimizing techniques, a 
trajectory entry time of 500 seconds was used for the thermal model . The convective 
heating rate history given in figure 2 is for a straight-wing shuttle vehicle entering 
at an angle of attack of 60°. The heating rate shown is for a 0.30S-meter-radius 
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(1- foot- radius) sphere. This heating rate was multiplied in the numerical program 
by a heating factor of 0 . 41 to provide surface temperature typical of the leading edge 
and fuselage nose areas . Although this is no longer the design ;ntry trajectory for 
the Space Shuttle , it is adequate for the comparisons required in this study . 


RESULTS 


The most efficient minimizing technique for all three TPS models was Peckham's 
method of least squares. For configuration 1 (no airgap). 12 function evaluations 
were required for convergence . For both configuration 2 (airgap between the Dyna- 
flex and the TG-15 000) and configuration 3 (airgap between the TG-15 000 and the 
aluminum) , 11 function evaluations were required for convergence. This represents 
a reduction in computer time by a factor from 13 to 16 when compared to the adaptive 
creep-pattern search method . The results for all cases are summarized in table 1 . 

The results for the adaptive creep-pattern search (ref. 1) are included in table I to 
assist in comparing the relative efficiency of each method . 

An ablative TPS model was also considered to determine whether P» ckham's 
technique would also be effective for ablators. The ablative TPS compot.ite model 
consisted of four materials: the Apollo ablator (Avco 5026/39/HCG), TG-15 000, 

titanium, and aluminum. The ablative TPS model was constructed as follows: 
ablator/titanium /TG-15 000/airgap/ aluminum. The control variables for optimization 
were the thicknesses of the ablator and the TG-15 000. For this case, the thickness 
of the titanium was 0.02 centimeter (0.008 inch) and the thickness of the aluminum 
was 0.18 centimeter (0.07 inch) . The maximum temperature at the backwall of the 
ablator was constrained to 922 ± 1 . 1 K (1660° ± 2° R) , whereas the backwall of the 
TG-15 000 was constrained to 422 ± 1 .11 K (760° ± 2° R) . No temperature constraint 
was imposed on the other two materials. The convective heating rate used was the 
same as for the passive TPS given in figure 1. This rate was modified by a 
factor of 0.8 to provide surface temperatures typical of the area of attachment of the 
wing lo the fuselage. The emissivity for the airgap was held at 0.8 on both the front 
and back of the materials . (The analytical ablation model and thermal properties are 
analyzed in reference 5.) For this case. 10 function evaluations were required for 
convergence, indicating that the efficiency of the least-squares routine for an 
ablating system is about the same as for a passive TPS system . 

The efficiency of using the quadratic fit for solving the one material and one 
constraint problem for both passive and ablative TPS problems was analyzed. In 
general , for a passive system , the solution was obtained with four function evalua- 
tions, whereas six or seven function evaluations were required for the ablative 
system. Because three function evaluations are required to provide the initial data 
for determining the coefficients of the quadratic ' the fourth evaluation for the 
passive system is used to confirm that the pre. oted thickness does produce the 
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desired temperature . The additional evaluations required for the ablative TPS are 
due to the nonlinearity induced by mass loss and its effect on the temperature distri- 
bution through the material . 


CONCLUSIONS 


The quadratic fit was the most efficient method for solving the one material and 
one constraint minimum weight problem for a thermal protection system . For the gen- 
eral problem involving more than one material and two or more constraints, Peckham's 
method of least squares was the most efficient minimizing technique investigated. 


Lyndon B . Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas. January 16, 1974 
986-15-31-04-72 
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TABLE I.- A COMPARISON OF THE RELATIVE EFFICIENCY OF DIFFERENT 
MINIMIZING TECHNIQUES FOR FINDING THE MINIMUM WEIGHT OF A 
THERMAL PROTECTION SYSTEM 


Method 

Number of function evaluations 

Configuration 1 

Configuration 2 

Configuration 3 

Adaptive creep^ 

191 

147 

171 

VA04A 

88 

79 

43 

VA02A 

44 

^29 

32 

MINSQ 

12 

11 

11 


Data are derived from reference 1 . 

Corrections became too small; converged in 15 more function evalu 
ations for a tot*>* of 45. 
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Figure 1 . - The three different TPS models used in the 
efficiency study. 
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